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Abstract 

We have investigated the recently proposed self-consistent theory of fluctuation - 
induced transport. In this framework the subsystem under study is coupled to two inde- 
pendent baths at different temperatures. In this non-equilibrium system one can extract 
energy at the expense of increased entropy. This is a simple model of Maxwell's demon 
engine that extracts work out of a non-equilibrium bath by rectifying internal fluctua- 
tions. We point out the errors in the earlier results. We have obtained an analytical 
expression for the fluctuation - induced transport current in a non-equilibrium state and 
various cases of physical interest have been elucidated. 
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From thermodynamics, it is well-known that useful work cannot be extracted from 
equilibrium fluctuations [1] . In a thermal equilibrium state the principle of detailed bal- 
ance ensures that no net particle current can flow in the presence of external potential 
of arbitrary shape. In contrast, in a non-equilibrium situation, where detailed balance 
is lost, net current flow is possible, i.e., one can extract energy at the expense of in- 
creased entropy. Several models have been proposed recently in this direction [2-6] . The 
motion of a heavily damped Brownian particle, in the presence of asymmetric static ex- 
ternal potential and under non- white or correlated fluctuations, is a simple example of a 
non-equilibrium system. In such a system induced current or directed motion appears, 
eventhough the average of the driving fluctuations vanishes. It turns out that the pre- 
ferred direction of motion and the magnitude of induced current depends sensitively on 
the shape of the potential as well as on the statistics of the fluctuations. Moreover, to 
obtain induced current the strength of non-white noise must exceed a certain minimum 
value and the magnitude of induced current shows a maximum value at an intermediate 
value of the noise strength. These models [2-7] of engines to obtain coherent response (or 
rectiflcation) from unbiased forcing come under the common denomination of "thermal 
ratchets" or "fluctuation - induced transport systems". The idea of thermal ratchets 
has been utilized recently for molecular separation [8]. One of the major motivations of 
these studies comes from molecular biophysics, where ratchet like mechanism is proposed 
to explain unidirectional movement of macromolecules or molecular motors. This is a 
physical example of preferred directional motion of Brownian particles (macromolecules) 
along periodic structures in the absence of obvious driving potentials, such as chemical 
potential gradients or thermal gradients. 

In a physically well motivated recent work [9] , Millonas points out that all the pro- 
posed earlier models [2-6] are basically phenomenological in nature and no attempt has 
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been made to formulate the problem from first principles. Millonas in his treatment 
[9] constructs a Maxwell's demon-type information engine that extracts work from non- 
equilibrium bath and allows a rigorous determination of kinetics consistent with the un- 
derlying laws of physics. He explicitly writes down a microscopic Hamiltonian including 
the subsystem and two thermal baths at different temperatures. An existing inequality 
of temperature can be exploited to do useful work. After eliminating bath variables one 
obtains non-linear Langevin equation for the subsystem variable Q : namely 

MQ + r(Q)Q + U'{Q) = U{t) + yJW)^B{t), (1) 

where T{Q) = {TA + ^Bf{Q))i iA{t) and iB{t) are two independent Gaussian white-noise 
fluctuating forces with statistics given by 

< U{t) >= 0, < U{t)U{t') >= 2 TAkT5{t - t'), (2.a) 

and 

< ^B{t) >= 0, < CB{t)CB{t') >= 2 TBkTSit - t'), (2.6) 

T and T are temperatures of the two baths A and B, respectively. Henceforth we set 
Boltzmann constant k to be unity. The bath B is characterised by a space dependent 
friction coefficient TsfiQ)- In eq.(l) U'{Q) is an external force (correspondingly U(Q) is 
an external potential) and we have used expression f(Q) for [V((5)]^ of ref.[9]. The final 
calculations in ref.[9] are done in an overdamped limit by simply neglecting Q term in 
eqn.(l). We have shown elsewhere that such a procedure is incorrect [10]. The resulting 
Fokker-Planck equation (eqn.(8) of ref.[9]), turns out to be inconsistent. For example, in 
an equilibrium situation (T — T), when the potential U(Q) is unbounded and positive, 
i.e., U{Q) — > oo as Q — >^ ±oo, the equilibrium distribution comes out to be 

P,{Q) = Cr(Q)e-^(Q)/^, (3) 
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where C is the normalization constant. This distribution function is incorrect because in 
equihbrium Pe{Q) must have the form 

PeiQ) = iVe-^W)/^, (4) 

where N is the position independent normalization constant. Also, if we set Fa — 
(i.e, the particle is coupled to a single bath B at temperature T) again the equilibrium 
distribution of the form of eqn.(3) is obtained. The reason for this inconsistency can 
be traced back to the improper overdamped limit of the original Langevin equation 
(eqn.(l)). For example in the absence of thermal bath A, Sancho et al. [11] have shown 
that the correct ovedamped Langevin equation should be 



rBf{Q)Q = -u\Q) - ^ 

The above equation leads to the correct equilibrium distribution as mentioned above. 
The error in the Millonas treatment [9] follows from the incorrect overdamped limit, 
where he ignores a term like the second term on the right hand side of eqn.(5). 

In this brief report, we construct the correct Fokker-Planck equation in the over- 
damped limit. We use this equation to study fluctuation - induced transport in a system 
where the function U(Q) and f(Q) are periodic under translation Q Q + X, {U{Q + X) = 
U{Q) and /(Q + A) — f{Q))- We obtain the correct expression for the mean velocity 
< Q > and study several special cases of physical interest. 

Following ref.[ll], after a straightforward algebra one can readily obtain Fokker- 
Planck equation for the variable Q(i.e, the evolution equation for the probability density 
P(Q,t)), in the overdamped limit and is given by 
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dt dQT{Q) ^ '^dQT{Q)dQT{Q) " 

""dQ T{Q) dQ T{Q) ^ ^dQ [T{Q)]^ ' ^ ^ 

This is the correct Fokker-Planck equation in the overdamped hmit and represents the 
Smoluchowski approximation to the original eqn.(l). When the potential U(Q) is un- 
bounded and positive, i.c, U{Q) ^ oo as Q ^ ±00, the system evolves towards the 
stationary distribution Pg{Q). This stationary distribution is characterized by no net 
current flow and is given by 

P,(Q) = iVe-*(Q), (7) 
where N is normalization constant and 

J \{TTA+TrBf(x)) r(x) {TTA+TTBfix))) ^' 

One can readily notice that in the equilibrium situation, i.e., when T = T, Ps{Q) reduces 
to the correct equilibrium distribution as given in eqn.(4). 

To study the case of fluctuation-induced transport, we take a simple case where both 
U(Q) and f(Q) are periodic functions and are invariant under the same transformation 
Q — > Q + A. Now, the basic problem reduces to flnding the mean velocity < Q > of the 
subsystem given the shape of U(Q) and f(Q). Following the procedure of refs. [9,12,13] 
closely, one can get the exact expression for the averaged velocity 

< Q >= [l-exp{-5)] 

where 

5 = *(a;) -*(a; + A), (10) 
5 



and ^ is given by eqn.(8). It is easy to see from eqns.(9) and (10), that in the equihbrium 
case, when the temperature difference between the baths is zero (T = T), the current 
vanishes indenticaUy (since 6 = 0). Also, one can easily varify that when the subsystem 
is coupled to a single bath, i.e., when either T a or Tb is zero, no net current is possible. 
It should be noted that the bath B which gives rise to space dependent friction coefficient 
TsfiQ) plays a special role. This can be noticed from the fact that if f(Q) is independent 
of Q the induced current is zero. In the extreme case of high friction limit (Ta ~^ 
CG or Vb oo) no net current is possible because the particle cannot execute the 
Brownion motion. 

There are two interesting special cases, one being the case of U(Q) = 0. It fol- 
lows from eqn.(9) that even in the absence of external potential fluctuation-induced 
currents are possible. In this case direction of the current depends on the details of 
the function f(Q). We see that the current will flow in one direction if T < T and in 
the opposite direction if T > T. However, their magnitudes are different. Thus, the 
system acts like a Carnot engine, which extracts work by making use of two thermal 
baths at different temperatures. In the other case where amplitude modulations in f(Q) 
and /'(Q) are small compared to the amplitude modulation of the potential U(Q), the 
second term in the eqn.(8) can be neglected. In this particular limit the problem be- 
comes equivalent to a particle moving in a spatially varying temperature - fleld, namely, 
T{Q) = (TrA + TrBfiQ))/T(Q). It is well-known in earlier literature [12-14] that such 
a spatial modulation of temperature-field can give induced currents. Moreover, the prob- 
lem of evolution of a particle in a spatially varying temperature field has a fundamental 
consequence in relation to local versus global stability criterion in statistical mechanics 
[15-16]. In a thermal equilibrium state the Boltzmann factor gives a relative occupa- 
tion probability of states at different local stable points without invoking the behaviour 
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of the potential profile between the states. However, in the presence of spatially vary- 
ing temperature- field, the relative stability (or occupation probability) between different 
states depend sensitively on the intervening potential and more importantly one can 
control the relative stability among two different state by modfying the kinetics in the 
sparsely occupied intervening state [14,16]. 

In conclusion, we have reinvestigated the recently proposed self-consistent micro- 
scopic theory of fluctuation-induced transport. We have pointed out errors in the earlier 
theory. We have obtained an analytical expression for fluctuation-induced current in a 
non-equilibrium situation and have discussed various cases of physical interest. Finally, 
we would like to mention that our expression for the stationary distribution (or steady 
state) Ps{Q) (eqn.(7)) is not a local function of U(Q) and f(Q). In such a situation as 
discussed above the relative stabiliy between two different local states in U(Q) depends 
sensitively on the intervening behaviour of U(Q) and f(Q). This can lead to interesting 
physics. For example, for given external potential U(Q), as one varies physical param- 
eters (in the parameter space of T, T, and Tb) we expect additional maxima and 
minima to appear in Ps{Q). Thus we can modify the stability properties of the sub- 
system. The qualitative changes in the stationary state of the subsystem is reflected in 
change in the behaviour of the extrema of the Ps{Q). Moreover, extrema in Ps{Q) may 
have little or no relationship to the extrema in the original potential U(Q). Each new 
structure in Ps{Q) may correspond to an entirely new state of a subsystem. In the spirit 
of the well known noise-induced phase transitions this is equivalent to having hierarchy 
of phase transitions, in a non-equilibrium system [17]. 
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